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Abstract 

in We perform a global analysis of the tree-level open-string amplitudes in the limit s <C Ma . 

o 

Based on the present data from the Tevatron, HERA, and LEP 2, we set a lower limit on the 

<N ■ 

£>y string scale Ms > 0.69 — 1.96 TeV at 95% confidence level for the Chan-Paton factors \T\ = — 4. 

a : 

We also estimate the expected sensitivities at the LHC, which can be as high as 19 TeV for \T\ = 4. 
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I. INTRODUCTION 

String theory provides a consistent framework for the unification of gravity with the 
standard model (SM) gauge interactions. The string scale Ms is naturally close to the 
quantum gravity scale M p i « 10 19 GeV. However, Arkani-Hamed, Dimopoulos, and Dvali 3] 
have proposed an alternative to this scenario. They assumed the space-time is more than 4 
dimensions with the SM particles living on a D3-brane. While the electromagnetic, strong, 
and weak forces are confined to this brane, gravity can propagate in the extra dimensions. 
The effective fundamental Planck scale is of the order of TeV with large extra dimensions 
as large as mm. TeV scale string theories are also possible with advance of D-branes jj, 
0). Thus, make it possible for experiments to test and probe the theory at high energy 
colliders [5|. 

In the simplest open-string model [(1 B, Isj] , it is assumed that all SM particles are identified 
as open strings confined to a D3-brane universe while a graviton is a closed string propagating 
freely in the bulk. One novel feature of the open-string model is the appearance of string 
resonances in the scattering of particles. Earlier works on phenomenological studies BBS 



lCj based on open-string amplitudes have been performed, which reduce to SM tree-level 
amplitudes at low energies . When s > Mj, clear resonances can be observed in experiments, 
in particular the mass spectrum of the resonances is given by \/nMs(n = 1,2,3...), which 



is very different from the Kaluza-Klein states in the ADD |2j or Randall- Sundrum [11 1 
model. On the other hand, when s < M| there would not appear resonances but may 
experience virtual interference effects from the string resonances. Effectively, the virtual 
effect is similarly to a contact interaction scaling as 1/M S . Note that the usual 4-fermion 
contact interaction due to exchange of new gauge bosons or leptoquarks has the leading 
effect of 1/A 2 , where A is the scale of new physics jl2j |. Bounds on four-fermion contact 
interactions induced by string resonances have been estimated based on the analysis on 
1/A 2 contact interaction jj. Although the approximation is valid, we shall show in this 
work that the effect is more complicated than just the simple substitution 1/A 2 — > 1/M S . 

The major improvement from previous work is as follows. We start directly from open- 
string amplitudes and per form a global analysis instead of using the result of an analysis 
on contact interaction |l3| . The global data include Drell-Yan production at the Tevatron, 
HERA charged and neutral current deep-inelastic scattering, total hadronic, and leptonic 
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e + e~ — > e + e~ , fi + fi~ , t + t~ cross section at LEP 2. Some of the data (CDF, HERA, and 
LEP 2) have been updated since the summer of 2004. We shall see that Drell-Yan production, 
due to large invariant mass data, provides the strongest constraint among the global data. 
By combining all data, the string scale Ms must be larger than 0.69 — 1.96 TeV for |T| = 0—4 
at 95% confidence level (CL). 

The organization of the paper is as follows. In the next section, we describe the tree-level 
open-string amplitudes and its low energy approximation. In Sec. IIHI we describe the high 
energy data sets that we used in this analysis. In Sec. IIVI we present our results on the fits 
and limits. In Sec. we estimate the sensitivity at the LHC We summarize in Sec. IVII 



II. OPEN-STRING AMPLITUDES 



tree-level open-string amplitudes follows the same approach as that 
They assumed that the tree-level open-string amplitudes 



Our construction o 
outlined in refs. 

represent the scattering of massless SM particles as the zeroth string modes. The SM 
is embedded in a type IIB string theory whose 10- dimensional space has six dimensions 
compactified on a torus with a common period 2ttR. There are M coincident D3-brane, on 
which open strings may end, that lie in the 4 extended dimensions. This approach assures 
the correct low energy phenomenology as given by the SM, yet captures one of the essential 
features of string theory, namely the string resonances. 

The 4-point tree-level open-string amplitude is given as a sum of ordered amplitudes 
multiplied by group theory Chan-Paton factors and can be expressed generically as 0, Q, 0] 

Astring = S(s, t)A 1234 T 12 34 + S(t, u)A 13 2 4 T 13 24 + S(u, s) A 1243 T 1243j (1) 

where A^^s are kinematic parts of the amplitude. The Mandelstam variables are denoted 
by s,t and u. Tiki's are the Chan-Paton factors Q| that involve traces over the group 
representation matrices A. For example, 

T1234 = MA1A2A3A4) +MA4A3A2AO. 

Typically, with normalization Tr(X a \ b ) = 5 ab the Chan-Paton factors are in the range of —4 
to 4. Since a complete string model construction for the electroweak interaction of the SM 
is unavailable, Chan-Paton factors are taken as free parameters in our calculation. S(s, t) is 
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essentially the Veneziano amplitude 



g(M H r V^-f> , P) 

r(l — a's — a't) 

where the Regge slope a' = M^ 2 , and S — > 1 as either s/M| or i/M| — ► 0. 

We start with Drell-Yan production of a pair of leptons. In the massless limit of the 
fermions we label their chirality a, (3 = L,R. The tree-level open-string amplitude is 



AstringiqLQL -> £r£r) = ig S 



Ti2uS(s,t)- +Ti 32 4S(t,u)- + Ti243S(u,s)(— - -) 

s u s u 



(3) 



The corresponding standard model amplitude is 



•Asm -> £r£r) = ig 2 -FRL , (4) 



where 



F a p = 2Q e Q q sin 2 6 W H ^"^f . (5) 

s — m z cos W 

Here = T 3/ - Q /S in 2 £ w , ^ = -Q /S in 2 £ w , Q f is the electric charge of the fermion / 
in units of proton charge, and the SU(2)i coupling g = e/sin# w . Identifying the string 
coupling with the gauge coupling g$ = g and matching the Chan-Paton factor as 



1243 — ^1324 



T, T 1234 = T + 1 (6) 



one can demand the string amplitude expression in Eq. (J2J) to reproduce the standard model 
amplitude in the low-energy limit. Since the string amplitude describes massless particle 
scattering, the masses of the W and Z gauge bosons are introduced by hand. So Eq. (J3J) 
becomes 

AstringiqLQj ^ £r!r) = ig z 'S(s,t)±F RL + ig 2 J '±-J \s ,t,u) , (7) 
f(s, t, u) = uS(s, t) + sS(t, u) + tS(u, s) . (8) 

The amplitude for qR~qk~ — > Il^l process is obtained from Eq. (jZJ by changing T R l to Tlr- 
As for the q R qR —> £r£r and qiJil — > £l£l processes, we use the corresponding T RR and 
J-~ll, respectively, and interchange t <-> u in Eq. (JJJ). 

The factor T(l — s/Ms 2 ) in the Veneziano amplitude S(s,t) develops poles at s — 
nM s {n = 1,2,3...), implying resonance states with masses y/nMs- In the limit Ms 3> 
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f(s,t,u) 



6 M| : 

7T 2 StU 



7T 2 



7T 2 US 



Therefore the amplitude squared for Drell-Yan process is given by 

Y, \M? = An 2 (\M e L l(s)\ 2 + \M"Us)\ 2 ) + 4t 2 (\M^ R (s)\ 2 + \M% L (s) 

where 



Mil, 



QeQq 



+ 



q e q q 

Hail a 



sin 2 9 W cos 2 9 W s - M| 



ix 2 uT 



7T 2 SU 



sin 2 Oy/Mg 



(9) 
(10) 



a = L,R 
(11) 



QeQq 



9a9} 



sin 2 6<w cos 2 9 W s - Ml 




ix 2 st 
6M| 



1T 2 tT 



4 sin 2 



a, j3 = L, R 
(12) 



Based on the above formulas the square of amplitudes for other processes can be obtained 
by crossing the Mandelstam variables and replacing the corresponding gauge boson. 

III. HIGH ENERGY PROCESSES AND DATA SETS 



Before describing the data sets used in our analysis, let us first specify certain important 
aspects of the analysis technique. Since the next-to-leading order calculations do not exist for 
the new interaction, we use leading order (LO) calculations for contributions from both the 
SM and the new interactions, for consistency. However, in many cases, e.g. in the analysis 
of precision electroweak parameters, it is important to use the best available calculations of 
their SM values, as in many cases data is sensitive to the next-to-leading and sometimes even 
to higher-order corrections. Therefore, we normalize our leading order calculations to either 
the best calculations available, or to the low-Q 2 region of the data set, where the contribution 
from the string resonances is expected to be vanishing. This is equivalent to introducing 
a Q 2 -dependent i^-factor and using the same i^-factor for both the SM contribution and 
the effects of string resonances. The details of this procedure for each data set are given in 
the corresponding section. Wherever parton distribution functions (PDFs) are needed, we 



use the CTEQ5L (leading order fit) set [16(. The reason to use the LO PDF set is that LO 
PDFs are extracted using LO cross section calculations, thus making them consistent with 
our approach. 



A. Drell-Yan Production at the Tevatron 

Both CDF 3] and D0 measured the differential cross section da/dM^ for Drell-Yan 
production, where Mm is the invariant mass of the lepton pair. (CDF analyzed data in both 
the electron and muon channels; D0 analyzed only the electron channel.) The differential 
cross section, including the contributions from the string resonances, is given by 

u y q 

where MJL is given by Eqs. (fTT|) and (|12p. s = M u , y/s is the center-of-mass energy in the pp 
collisions, Ma and y are the invariant mass and the rapidity of the lepton pair, respectively, 
and £1,2 = ^fe ±y . The variable y is integrated numerically to obtain the invariant mass 
spectrum. The QCD i^-factor is given by K = 1 + (1 + ^f~)- We scale this tree- level 
SM cross section by normalizing it to the Z-peak cross section measured with the data by 
a scale factor C (C is very close to 1 numerically). The cross section a used in the fitting 
procedure is given by 

a — C (a"SM + Cinterf + CT str i ng ) , (13) 

where cr irLterf is the interference term between the SM and the string resonance amplitudes and 
^string is the cross section due to the string resonances interactions only. When normalizing 
to the low-energy data, we neglect the possible contribution from the string resonances, as 
it is much smaller than the experimental uncertainty on the data that we use. 



B. HERA Neutral and Charged Current Data 



ZEUS 



2Cj and HI 2l( have published results on the neutral-current (NC) and charged- 



current (CC) deep-inelastic scattering (DIS) in e + p collisions at w 319 and 318 GeV, 
respectively. The data sets collected by HI and ZEUS correspond to an integrated luminosi- 
ties of 65.2 (HI), 60.9 (ZEUS CC) and 63.2 (ZEUS NC) pb" 1 . ZEUS 0,Q and HI 24 1 
have also published NC and CC analysis for the data collected in e~p collisions at y/s ^318 
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and 320 GeV, respectively, with an integrated luminosity of 16.4 (ZEUS CC), 15.9 (ZEUS 
NC), and 16.4 (HI) pb" 1 . 

We used the single differential cross section da/dQ 2 presented by ZEUS Q and 
double differential cross section d 2 a/dxdQ 2 published by HI 2J|. The double differential 
cross section for NC DIS in the e + p collisions, including the effect of the string resonances, 
is given by 

d 2 a 



(e + p -> e + X) 




+ 



5> 



x) 



vn\M? L (t)r + \M e A(m + i^wr + wzm 



mi(t)\ 2 + \MZ(t)\ 2 + (i - yn\Mi« R (t)\ 2 + wzm 



(14) 



where Q 2 = sxy is the square of the momentum transfer and f q /q(x) are parton distribution 
functions. The reduced amplitudes M*% are given by Eqs. (fTTj) and (|T2*|). The double 
differential cross section for CC DIS, including the effect of string resonances, can be written 

as 



d 2 a 
dxdQ'- 



(e + p^uX) 



9 

64vr 



7T til 

6Mf 



TT 2 St 





2 


\ n 2 uT 




) 2M| 






2 


\ 7T 2 ST 




/ 2M| 





(l-y) 2 (rf(x) + s(x)) 
+ c(x) 



(15) 



where <i(x), s(x),u(x), c{x) are the parton distribution functions. The single differential cross 
section da/dQ 2 is obtained from the above equations by integrating over x. The cross section 
in the e~p collisions can be obtained by interchanging [LL <-> LR, RR ^ RL) in Eq. ()14|) 
and by interchanging (q(x) ^(x)) in Eq. fl3j) . 

We normalize the tree-level SM cross section to that measured in the low-Q 2 (Q 2 < 2000 
GeV 2 ) region. The cross section used in the fitting procedure is then obtained similarly to 
that in Eq. (fHl) . 



LEP 2 Data 



We analyze the LEP 2 observables sensitive to the effects of the string resonances, includ- 
ing hadronic and leptonic cross sections. The LEP Electroweak Working Group combined 
the qq,fi + fi~, and t + t~ data from all four LEP collaborations 2s| for the machine energies 



between 130 and 207 GeV. We use the following quantities in our analysis: (i) total hadronic 
cross sections, (ii) total t + t~ cross sections. We take into account the correlations 

of the data points in each data set as given by 2jjj. For Bhabha scattering cross section 
a(e + e~ —>■ e + e~), we use various data sets from individual experiments 
The angular distribution for e~e + — > f f (/ = q, e, fi, r) is given by 

da 



26 



dcos9 
N f s 

1287T 



1 + COS I 



M e / L (s) 



+ W R f R {s) 



COS I 



M e L f R (s) 



\M e J L (s) 



+5, 



<7 



costf) 2 (\M&{ 8 ) + Ml{{t)\> + \M^ R {s) + M^(t)| 2 - \Ml{{s) 



\M e R f R (s) 



+4(|M^(t)| 2 + |M^(t)| 



ef 



where Nf = 1 (3) for i (q), and M e L is given by Eqs. (jlljl and (|12j). The additional terms 
for f = e arise from the t,w-channel exchange diagrams. 

To minimize the uncertainties from higher-order corrections, we normalize the tree-level 
SM calculations to the NLO cross section, quoted in the corresponding experimental pa- 
pers. We then scale our tree-level results, including contributions from the Z, 7, and string 
resonances with this normalization factor, similar to Eq. (|13|) . When fitting angular dis- 
tribution, we fit to the shape only, and treat the normalization as a free parameter of the 
fit. 



IV. CONSTRAINTS FROM HIGH ENERGY EXPERIMENTS 

In the previous section, we have described the data sets from various high energy ex- 
periments used in our analysis. Based on the above individual and combined data sets, we 
perform a fit to the sum of the SM prediction and the contribution of the string resonances, 
normalizing our tree-level cross section to the best available higher-order calculations, as ex- 
plained above. As seen from Eqs. (jll)) and (|12jl. the effects of the string resonances always 
enter the equations in the form 1/M|. Therefore, we parameterize these effects with a single 
fit parameter 7]: 

1 
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In most cases, the differential cross sections in presence of the string resonances are bilinear 
in 77. We use MINUIT to minimize the \ 2 - 

The best-fit values of 77 for each individual data set and their combinations are shown in 
Table 01 In all cases, the preferred values from the fit are consistent with zero, and therefore 
we proceed with setting limits on 77. The one-sided 95% CL upper limit on 77 is defined as: 

0-95 = J ° , p , " , 16 
Jo d vP{V) 

where P(rj) is the fit likelihood function given by 

P( V ) = -±—ex P (-( X 2 ( V ) - J&J/2) . (17) 

The corresponding upper 95% CL limits on 77 and lower 95% CL limits on Ms are also shown 
in Table|U The combined limit is as high as 1.96 TeV (1.92) TeV for T = 4(— 4). Note that 
even for T = 0, there is still some stringy effects (1/Mj), as shown in Eqs. ([lip and (fT2*j) . 
It is clear from the results of Table H] that the dominant set of data affecting the limit on 
Ms is the Drell-Yan data. This is simply because the effect of string resonances scales like 
1/M|. We show the effect of string resonances on the invariant mass spectrum in Fig. ^ 

V. SENSITIVITY AT THE LHC 

We expect that at the LHC the best channel to probe the string resonances is Drell-Yan 
production. We assume that the approximation Mg 3> s, \t\, \u\ is still valid at the LHC. 
Therefore, the reduced amplitudes of Eqs. (fTTj) and (|T2j) can be used and tested by a direct 
comparison at the LHC. 

It was shown in Ref . [28^ that using the double differential distribution d 2 cr/ Mud cos 6 
can increase the sensitivity to the KK states of the graviton compared to the use of single 
differential distributions. Similarly, we expect this to be the case for the string resonances. 
The double differential cross section for Drell-Yan production, including the interactions of 
the 7, Z, and string resonances, is given by 



dM u dyd cos 9* ^ 192ns 



:i+coser(\Mii(s)\ 2 +\M e us)\ 2 ) 



:i-cos6r{\Mll(s)\ 2 + \M% L (s)\ 2 ) 
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100 200 300 400 500 600 700 800 900 
M ee (GeV) 

FIG. 1: The differential distribution da/dM ee for Drell-Yan production at the Tevatron (Run I) 
for the low scale open-string model and for the SM. 



where are given by Eqs. (fTT|) and (fT2")l . 9* is the scattering angle in the center-of-mass 
frame of the initial partons, s = M^, dx\dx2 = (2M«/ s)dMudy, and Xi % = Mae ±y / \/s. 

We briefly describe the procedures here, (i) We calculate the SM cross section in each 
bin. (ii) In each bin, we obtain the expected number of events by multiplying the SM cross 
section by the integrated luminosity, (iii) In each bin, we generate the number of events 
according to the expect number of events. If the expected number of events is less than 100, 
we used the Poisson statistics otherwise we used the Gaussian, (iv) After generating the 
data set, we fit the data set to rj using Eq. f)13|) . Then we evaluate the 95% CL lower limit 
on Ms implied by the data, (v) We repeat the above procedure for 1000 times. Then we 
histogram the limit obtaining each time. We take the sensitivity to be the medium (501st) 
of the histogram. 

In each simulation, we always normalize to the Z-peak data, by which a lot of systematic 
uncertainties are eliminated. We have used an efficiency of 90% and a rapidity coverage of 
\rj\ < 2.0 for lepton detection. We also used a single experiment with combined e, fi samples. 
The sensitivity may be lowered by about 10% if some level of systematic uncertainties are 
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used in calculating the \ 2 |2&|- The sensitivity, at the 95% CL, to M s at the LHC (100 fbr 1 ) 
with different Chan-Paton factors is given in Table |H] The sensitivity reach is as high as 19 
TeV for \T\ =4. Note that the sensitivity obtained in Ref. is somewhat lower than the 
sensitivity obtained in this paper, because their analysis is based on the search for a single 
resonance state while ours is based on the deviation in the whole invariant mass spectrum. 
Therefore, statistically we can obtain a higher sensitivity. 



VI. CONCLUSIONS 



The ultimate consequence of a low string scale is the testable string scattering amplitudes 
at hadron collider. The scattering of particles shows the stringy nature. The most obvious 
sign is the presence of string resonances. Even if the energy is not high enough, some 
deviations from the SM due to the string resonances are still expected. For example, the 
tree-level open-string amplitude for dilepton production at low energy can be expressed 
generically as 

Astring ~ A S Ai{s,t,u) ■ S(s,t,u) + Tf(s,t,u) ■ g(s, t,u) , 

where Asm is the SM amplitudes, S(s,t,u) are the Veneziano amplitudes given in Eq. ©, 
T is the undetermined Chan-Paton parameter, f(s, t, u) a kinematic function given in Eq. 
(fTUj) . and g(s,t,u) some process-dependent kinematic functions. In the low energy limit 
s <C Mj ,S(s,t) —> 1, f(s,t,u) — > 0, the open-string amplitude reduces to the standard 
model amplitude. 

In this work we have performed a global analysis of the model based on the data from 
Drell-Yan production at the Tevatron, HERA neutral and charged current deep-inelastic 
scattering, and fermion pair production at LEP 2. Drell-Yan production at the Tevatron, 
due to large invariant mass data, provides the strongest constraint. Combining all data the 
string scale M s must be > 0.7 - 2 TeV for Chan-Paton factors \T\ = - 4 at 95% CL. We 
have also estimated the expected sensitivity at the LHC, which is about 19 TeV for \T\ = 4. 
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TABLE I: Best-fit values of rj = 1/Mj and the 95% CL upper limits on r] for individual data set 
and combinations. Corresponding 95% CL lower limits on Ms are also shown. By default we use 
T = 4 in part(a). We show in (b) limits for all combined data with different T. 
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TABLE II: Sensitivity to the parameter r/ = 1/M S at the LHC, using the dilepton channels. The 
corresponding 95% CL lower limits on Ms are also shown. 



LHC (14 TeV, 100 fb" 1 ) Dilepton 

T=4 

T=l 

T=0 
T=-l 
T=-4 



^5 (TeV 4 ) 
8.49 x 10~ 6 
3.71 x 10~ 5 
5.79 x 10~ 4 
3.42 x 10~ 5 
7.95 x 10~ 6 



95% CL lower limit on M s (TeV) 
18.5 
12.8 
6.4 
13.0 
18.8 
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